General framework for estimating the ultimate precision limit in noisy quantum-enhanced 

metrology 
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The estimation of parameters characterizing dynamical processes is central to science and technology. The 
estimation error decreases with the number N of resources employed in the experiment (which could quantify, 
for instance, the number of probes or the probing energy). Typically, it scales as 1/s/N. Quantum strategies 
may improve the precision, for noiseless processes, by an extra factor l/\f~N. For noisy processes, it is not 
known in general if and when this improvement can be achieved. Here we propose a general framework for 
obtaining attainable and useful lower bounds for the ultimate limit of precision in noisy systems. We apply this 
bound to lossy optical interferometry and atomic spectroscopy in the presence of dephasing, showing that it 
captures the main features of the transition from the 1 /N to the 1 / yN behavior as N increases, independently 
of the initial state of the probes, and even with use of adaptive feedback. 



It is by now well established that, in the absence of noise, 
quantum properties of the probes help to increase the pre- 
cision in the estimation of parameters that characterize dy- 
namical processes, like phase shifts in optical interferometry 
or transition frequencies in atomic spectroscopy fl-Q]. The 
estimation error decreases with the amount of resources em- 
ployed in the measurement, which might be for instance the 
energy of a probing light field or the number N of identical 
probes. For independent probes, it is proportional to 1/yN, 
a consequence of the central-limit theorem. By entangling the 
probes one may attain, for noiseless processes, the ultimate 
lower bound for the estimation error, which scales then with 
l/iVjSlilh the so-called Heisenberg limit, thus allowing bet- 
ter accuracy for the same number of resources. For noisy pro- 
cesses, however, it is not known in general if this limit can be 
attained, and if entanglement can still be a helpful resource for 
this purpose. For lossy optical interferometry, it was shown 
recently |5|,|6|] that the Heisenberg limit is not attainable when 
N — » oo. General expressions for the uncertainty in the esti- 
mation are known, but their calculation involves complex opti- 
mization procedures, which become quite cumbersome when 
the number of resources increases. Here we show that the ef- 
fectiveness of quantum states for parameter estimation in the 
presence of noise can be precisely assessed. We introduce 
a bound for this uncertainty, proven to be attainable, which 
leads to useful expressions for the ultimate limit of precision 
in noisy systems. We exemplify the utility of this bound by 
applying it to lossy optical interferometry and atomic spec- 
troscopy in the presence of dephasing. 

Parameter estimation and the Cramer-Rao bound 

A general protocol to estimate an unknown parameter x cor- 
responding to a quantum process is shown in Fig. la). This 
parameter is estimated from the knowledge of the initial and 
final states of a probe that undergoes the process under investi- 
gation. The protocol may be divided in three stages. First, the 
probe is prepared in an initial state (first box on the left-hand 
side of the figure), and evolves under the action of the quan- 
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FIG. 1: Setups for quantum parameter estimation, a) General 
algorithm to estimate an unknown parameter x of an arbitrary dy- 
namical process. The probe, prepared in a known initial state, is sent 
through a physical channel. A measurement is performed on the final 
state, from which the parameter x is estimated, b) Setup for estimat- 
ing a phase shift 6 in an optical interferometer. 



turn process, corresponding to the second box in Fig. la). In 
general, after this step, the probe will be in a mixed state (third 
box). The second stage concerns the choice of a suitable mea- 
surement, applied to the probe after its evolution (fourth box), 
so as to extract information about the parameter to be esti- 
mated. The third and last step consists in associating, through 
some rule (estimator), each experimental result with an es- 
timation of the parameter. A concrete example is shown in 
Fig. lb), which refers to phase shift estimation in optical in- 
terferometry. 

Typically, one has some prior knowledge about the param- 
eter. This is sometimes called in the literature "local estima- 
tion", as opposed to the situation when there is complete ig- 
norance about the parameter to be estimated, which is known 
as "global estimation." In local estimation, a merit quantifier 
that allows the comparison of different strategies and perfor- 
mances is the error estimate 10, d] fix = J ((x cs t — £reai) 2 ), 
where the average is taken over all possible experimental re- 
sults, x rca i is the real value of the unknown parameter x and 
x es t is the estimated value of x, obtained from the measure- 
ment results through the use of an estimator. 

When the initial state, the physical system, and the mea- 
surement procedure, involving v repetitions of the experiment, 
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are fixed, the error 5x is limited by the Cramer - Rao inequal- 
ity |0] Sx > 1/ y / vF(x rea i) (valid for unbiased estimators, 
that is, estimators for which (x cst ) = x rca i). Here F(x) is the 
Fisher Information, given by 



d\n[ Pj (x)} 
dx 



(1) 



where pj (x) is the probability of getting the experimental re- 
sult j. This relation is valid for both classical and quantum 
physics. 



In quantum mechanics, pj (x) = Tr p(x)E 



where the 



operators Ej correspond to a specific measurement setup, as- 
sociated to the results j, and p(x) is the density matrix of the 
probe after its interaction with the system under investigation. 
The Hermitian operators Ej are positive operator-valued mea- 
sures (POVM's) [10], satisfying the relation £\ Ej = 1. In 
general, the evolution of a density matrix p under the action of 
a quantum channel may be expressed in terms of Kraus opera- 
tors ifTHl flg(x) as p(x) — J2e n^(x)poftJ(x), where p is the 
initial value of p and ^^n](x)n^(x) = 1, The maximiza- 
tion of (fl]i over all possible measurement strategies yields the 
so-called Quantum Fisher Information Fq SSI El and 
the quantum generalization of the Cramer-Rao inequality, 



SX > 1/^JvFq [p (x rC al)] 



(2) 



For a closed system, evolving under a unitary transforma- 
tion U (x) and prepared in a pure initial state po = \ip)(ip\, Fq 
can be expressed as IU411 



U(x)p U\x) 



= A(AH 2 



(3) 



where 



(AH 2 



(V|tf 2 (zM - ^\H{x)\^f 



and H (x) = i \dfj^{x)/dxj U(x). However for more gen- 
eral situations where the initial state is a mixture and/or the 
evolution is not governed by a unitary transformation, a closed 
expression for Fq as a function of Kraus operators is un- 
known and finding it remains an open problem. 

Although analytic expressions for Fq have been found for 
some specific classes of initial states and non-unitary pro- 
cesses (see, for instance refs. |[T5l - H7ll ). in most cases only 
an upper bound to Fq can be obtained. For example, when 
the initial state of the probe is mixed but the quantum chan- 
nel is unitary, an upper bound based on the convexity of Fq 
can always be established. More recently an upper bound to 
Fq, applicable to general quantum channels, was derived in 
terms of Kraus operators and was shown to be attainable for a 
special class of processes ("quasi-classical" processes) fljll . 
Bounds for error estimation in noisy systems 
We introduce here an universal upper bound to Fq, expressed 
in terms of the initial state of the probe and any Kraus repre- 
sentation of the quantum channel. This bound is valid for both 



unitary and non-unitary processes, and it is always possible to 
choose a Kraus representation such that it coincides with the 
quantum Fisher information. We show that this bound leads to 
useful analytical approximations for the precision of parame- 
ter estimation. 

Let S be the probe used to estimate a parameter of a gen- 
eral dynamical process. Our strategy consists in introducing 
additional degrees of freedom, which play the role of an en- 
vironment E for the system S, so that a general dynamical 
evolution concerning S is transformed into a unitary evolu- 
tion for S + E. The problem is thus reduced to parameter 
estimation for a unitary evolution. A given general dynami- 
cal evolution of a system S may be related to an infinitude of 
unitary evolutions of an enlarged system consisting of S plus 
some "environment" E. Fixing the unitary evolution and the 
"environment" E is equivalent to singling a specific represen- 
tation of Kraus operators out from the unlimited number of 
such representations that describe the same dynamical evolu- 
tion of S alone. 

The detailed calculation of this upper bound can be found 
in the Supplementary Material. An outline of the derivation is 
presented in the Methods section. Given an initial pure state 
po = \ip)(ip\ of a system S and an arbitrary process, which 
changes the state to p(x) = J2e ri^(x)/5ori](x), we get the 
upper bound to Fq 



where 



C Q (po,Ue(x))=A (ffi(x)) - (H 2 (x)) 



Hi{x) = V ( x ) jft* ( x ) 

' dx dx 



H 2 (x) 



dfl\ (x) 



dx 



Il e (x) . 



(4) 



(5) 



(6) 



the symbol ( ) meaning (•) = Tr [•po]. 

Equations (HJ)-© show that the upper bound Cq can be ex- 
plicitly evaluated in terms of the Kraus operators describing 
the quantum channel and the initial state of the probe. When 
the process is unitary Cq reduces to the quantum Fisher in- 
formation, as given by Eq. (01. This upper bound can also be 
applied to situations with ancillas, where the initial state of the 
probe is entangled with some external system. The bound is 
then easily calculated from the reduced density matrix of the 
probe. The attainability of this bound is demonstrated in the 
Methods section, implying that 



■Fq [p{x)\ = mm {fle(x)} C Q p ,Iie(x) 



(7) 



where the minimization runs over all Kraus representations 
{Il£(x)} of the quantum channel. 

Equation (|7]) leads to an interesting physical insight, within 
the framework of our method: there is always an environment 
E such that monitoring it together with the system S does 
not lead to more information about the parameter x than that 
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obtained by monitoring just the system S itself. These con- 
siderations motivate a strategy for choosing convenient Kraus 
operators: the aim should be to reduce the non-redundant in- 
formation about the parameter in the environment. 

We apply now the bound dU to some important problems 
in quantum metrology, and show that it leads to the elucida- 
tion of fundamental questions regarding limits of precision in 
parameter estimation. 

Role of entanglement in quantum-enhanced metrology 

In recent years, much effort has been made to establish 
the relevance of entanglement in quantum metrology |1- 
lil [l6l fnl [l9l - l25ll . For unitary processes, it is now clear that 
entanglement does lead to enhancement of the precision in pa- 
rameter estimation. However, an important question remains: 
does this gain persist in the presence of decoherence? 

A typical paradigm of quantum metrology involves N 
probes undergoing N independent and identical x-dependent 
processes, and being submitted to a measurement afterwards, 
with the aim of estimating the parameter x. If the initial state 
of the probes is not entangled and each probe is measured 
independently of the others, then the error scales at most as 
Sx ~ 1/VN h. It was shown in [4] that, for unitary pro- 
cesses and in the absence of feedback, initial entanglement is 
a necessary ingredient to improve the error scale to the ulti- 
mate quantum limit Sx ~ 1 /N. Moreover, it was also shown 

that, for unitary processes of the type U(x) = exp (ixH^J, 

it is always possible to find an initial entangled state of the 
probes that leads to Sx ~ 1/N. 

We consider now the effect of noise on these results. When 
the N probes are submitted to N independent processes, 
which depend on the parameter x, a set of Kraus operators 
that describes the evolution of the probes is 



n flA ,...,^(x) = ir^(x) ® u^'(x) ® ... ® u^>(x), (8) 

where Tig (x) is a Kraus operator corresponding to the m-th 
process. From (|4j, it is then straightforward to show that Cq 
can be decomposed into two parts (see Supplementary Mate- 
rial) 



rr(2), 



JY 



Cq (/> : ,H, : ,, :.r:) = 4 E [<^ m) > - (H ( 2 m) ) 2 

m=l 

N mi-1 

+ 8 E E [(^ mi) ^ m2) )-(^ mi) x^ m2) ) 

7711— 2 77*2=1 



,(9) 



where the operators h["^ correspond to the definitions Q 
and © for the m-th process. The first term in Eq. (O is of 
o(N), while the second is of o(N 2 ), as they involve the sum 
of N and N(N — l)/2 terms respectively. 

Since Cq > Fq, this result has two immediate implica- 
tions: (i) a necessary condition for Sx ~ 1/N is that the sec- 
ond line of (O is different from zero, which establishes the 
kind of correlation needed in the initial state of the N probes 
in order that the Heisenberg limit is attained - this condition 
is also sufficient if there is a Kraus representation in the form 
dHJ that minimizes Cq; (ii) if a quantum channel has some 



Kraus representation for which = 0, then Sx ~ 1/ VN 

at most, even when the initial state is entangled. The fact that 

H. 



implies Sx ~ 1/ v N was derived for the more 



(to) 

2 

restrict situation of finite-dimensional spaces in Il26ll . which 

also showed that the condition = holds for almost 

all finite-dimensional quantum channels, the "full-rank" chan- 
nels. Moreover, (ii) remains valid when feedback is included, 
as demonstrated in the Supplementary Material. Therefore, 
for this class of channels, the limit Sx ~ 1/N cannot be 
achieved, in the presence of decoherence, even with the use 
of entanglement and/or feedback control. 
Precision limits for lossy optical interferometry 
Optical interferometry with standard light sources leads to an 
uncertainty in the determination of the phase that scales with 
the inverse of the square root of the mean number of photons, 
the so-called shot-noise or standard quantum limit 1271 12811 . 
On the other hand, it has been shown that squeezed 11271, 12811 
or entangled states yl Hfl may lead, in the absence of losses, to 
a scaling of this uncertainty with the inverse of the mean num- 
ber of photons. This is the ultimate limit imposed by quantum 
mechanics II41 I29I1 . 

In the presence of losses, general results concerning the 
limit of precision are still unknown. Special entangled quan- 
tum states have been proposed for optical metrology in order 
to mitigate the deleterious effect of noise on the precision of 
phase estimation [23]. Numerical limits for the minimum un- 
certainty in the measurement of the phase shift, when either 
one or both arms of an interferometer are subject to photon 
losses, were found in il^[l7ll for input states with fixed pho- 
ton numbers. These studies show that, for states with total 
number of photons up to N = 80, the best possible scaling 
of the uncertainty is intermediate between the standard and 
the Heisenberg limit. Experimental confirmation with two- 
photon entangled states was presented in J[25fl. while bounds 
found in [5] for global estimation and in [6] for local estima- 
tion lead to a 1/yN scaling. This raises an important ques- 
tion concerning the behavior of the precision as a function of 
iV.Here we answer this question in a general way, by finding 
an analytical lower bound for the phase uncertainty. 

We consider for definiteness a two-arm interferometer with 
a dispersive object placed in the upper arm (see Fig. QJb)- An 
incoming photon is described by a two-mode state, each mode 
corresponding to one of the arms of the interferometer. In 
the absence of losses the initial two-mode state in the inter- 
ferometer evolves into \ip(9)) = U(0)\ip)o — exp(in9)\ip)o, 
where h is the photon number operator corresponding to the 
dispersive-arm mode and 9 is the phase shift parameter to be 
estimated. Therefore, from ([3), 



Fq (p ) = 4(An 2 ) and 89 > (4^(An 2 



-1/2 



, (10) 



where (An 2 )o is the photon-number variance of the state 
in the dispersive arm of the interferometer, before it under- 
goes the phase shift. From this expression, it follows that 
states maximizing the photon-number variance lead to min- 
imum values of the phase uncertainty. This is the case of 
the NOON states |^(A0) = (1-^,0) + \0,N))/y/2 UM, 
where |iVi , AT 2 ) represents a Fock state with Ni and N2 pho- 



4 




500 1000 



N 



FIG. 2: Lower bounds for the phase error. Lower bound 88 of the 
normalized variance 80^^l£- as a function of TV for 77 = 0.5 (blue), 
77 = 0.9 (green), r\ — 0.99 (red), and r\ = 1 (black). The respective 
dashed lines correspond to the limit d!2t for each value of 77. 




FIG. 3: Numerical test of the tightness of the bound. Comparison 
between the numerical maximum value of Tq and the upper bound 

Cq, \JFqICq as a function of 77, for TV = 10 (blue), TV = 20 
(green), and TV = 40 (red). The inset displays the minimum of this 
ratio over all values of 77 as a function of TV. 



tons in the upper and lower arms, respectively. For this state, 
(An 2 ) N ooN = TV 2 /4, which yields the Heisenberg limit 
86>\/{y/VN). 

We consider here interferometers with losses only in the 
dispersive arm - see Fig. lb). Generalization to losses in both 
arms is presented in the Supplementary Material. As shown 
in the Methods section, proper choice of Kraus operators, in- 
spired by physical considerations, leads to 



4?7(n)o(An 2 )o 



(An 2 ) (1 - 77) + r/(n) 



(11) 



where 77 quantifies the photon losses (from 77 = 1, lossless 
case, to 77 = 0, complete absorption) and (h)o is the initial 
average number of photons in the dispersive arm. 

The above equation has a high degree of generality: it 
was derived under no assumption whatsoever about the ini- 
tial state, which might therefore have a definite, limited or un- 
limited number of photons. Furthermore, it displays the limit 
of precision as a continuous function of losses as well as the 
variance and the average photon number of the field in the 
dispersive arm, for the entire range of values of these quanti- 
ties. When (Ah 2 )o/ (n)o << 77/(1 — 77), this expression ap- 
proaches the lossless case (fTOl.while when (An 2 )o/ (n)o >> 
77/(1 — 77) Eq. ( fTTT i yields the more restricted result 



> 



/ 1 -77 
4z/r?(n.) ( 



(12) 



Equations ( fTTb and dT2b imply, quite generally, that photon 
losses in interferometers gradually blur the gain yielded by 
the special quantum states for phase estimation; even with the 
best strategy, asymptotically the improvement with respect to 
standard light sources is not by a scale change but only by a 
limited constant factor, as shown in ( fl~2b . This last expression 
coincides with the bounds obtained in Q5J, |^] for states with 
definite photon number. 



For states with well-defined number TV of photons, exten- 
sively discussed in quantum metrology, Eq. (fTTT l leads to an 
upper bound for the quantum Fisher information maximized 
over all initial states, which depends only on 77 and on TV. In 
this case, one gets then (see Supplementary Material) 



J-Q(po) < C Q = 



2TV 



(i-v)N 

V 



(13) 



The behavior of the corresponding lower bound for the phase 
uncertainty 89, as a function of TV, is shown in Fig. [3] where 
the phase estimation error is normalized by the shot-noise 
limit, expressed in terms of the total number of photons in 
both arms Ql El: 89 SN = (y/fj + l)/y/4wqN. This re- 
sult clearly exhibits the change from the Heisenberg scale to 
a 1/vTV scale. For TV » 77/(1 — 77), one gets a 
scaling, corresponding to Eq. (fl2l l. with the average number 
of photons replaced by TV, while for TV << 77/(1 — 77), the 
Heisenberg dependence, proportional to 1/TV, is obtained. 

One may wonder how tight is the bound ( fT3T >. The com- 
parison of this bound with Tq for the numerically determined 
optimal states is shown in Fig. [3] for all 77 and for different 
values of TV up to TV = 100. This yields, for this range of TV, 



1/ Y vCq < 89 < 1.25/y vCq, thus showing that, for these 

states, Cq provides a very good qualitative and quantitative 
approximation to the ultimate quantum limit. As N increases, 
the fundamental limit of the phase uncertainty given by the 
quantum Fisher information becomes at most 1 .25 larger than 
the one given by our bound. 

Precision limits for atomic spectroscopy under dephasing 

We consider now the estimation of transition frequencies in 
atomic spectroscopy, in the presence of Markovian dephasing, 
which is a common source of decoherence for atoms, and the 
most important one in trapped-ions experiments. The aim is to 
estimate the transition frequency ujq, by preparing the atoms 
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in a known initial state, letting them evolve freely, and then 
measuring the final atomic state. The resources here, which 
affect the error in the estimation, are the number N of atoms 
and the total time T = tu, where t is the evolution time for 
each atom and v is the number of experimental repetitions. In 
the absence of decoherence, the error scales as ~ 1 j\f~N if the 
initial state of the N atoms is separable. On the other hand, an 
initial maximally entangled GHZ atomic state may improve 
the scaling to ~ 1/N IW. The situation changes when de- 
phasing is present. In 11911 . it was shown that, in this case, a 
separable state or a GHZ state of the atoms lead to the same 
error Sujq = ^J2^e/ (NT), where 7 is the dephasing rate. It 
was also shown that a generalized Ramsey spectroscopy, in- 
volving a joint measurement of all the atoms, may lead, for an 
optimal non-maximally entangled state, to an error limited by 
5ujq > v/27/ (NT), a better bound then before, but still pro- 
portional to 1/ y/N. In this regard, two questions were raised 
in Ref. [19]: whether this bound is saturated asymptotically 
by the proposed measurement strategy; and whether it coin- 
cides with the one obtained from the quantum Fisher infor- 
mation. Ref. IBOIl has shown that, for properly chosen initial 
states, this bound is actually saturated asymptotically, for the 
proposed measurement scheme. Here we show that the re- 
maining question can be tackled by our method. Indeed, we 
show that the above bound coincides asymptotically with the 
one given by the quantum Fisher information optimized over 
all initial atomic states. 

We show in the Supplementary Material that it is possible 
to choose a physically-motivated set of Kraus operators such 
that our bound, optimized over all initial states, leads to 



< 



NT 
17 



2jtN 



l+(e 



2jt 



l)N 



< 



NT 
^7 



(14) 



which coincides with the bound derived in Ref. [19]. Since 
Ref. 13011 showed that this bound is attained asymptoti- 
cally for the above generalized Ramsey measurement and 
a specific class of initial states, it follows that necessarily 
limjv^oo vFq &x /N > T/2-f. In view of inequality ([T4]i, one 
concludes that lirriAr^oo vTq^/N = T/2j. This general- 
izes the previous results, showing that, no matter the initial 
state and the measurement scheme, the maximum improve- 
ment obtainable in the presence of dephasing is by a factor 
^fe, which does not change the scaling of the error with N. 
Therefore, also in this case decoherence does not allow attain- 
ing the quantum scale, no matter how small the dephasing rate 
is, as long as N is sufficiently large. It can be shown that this 
conclusion remains valid in the presence of feedback. 
Summary and perspectives 

We presented here a lower bound for the error in single- 
parameter estimation, within the framework of quantum 
metrology, valid for both unitary and non-unitary processes, 
which is always attainable. The calculation of the best state 
leading to the ultimate quantum limit in the presence of noise 
is a difficult task, requiring numerical analysis that becomes 
more and more cumbersome as the number of resources in- 
creases. The upper bound |@), on the other hand, circumvents 
this difficulty, leading to useful relations that do not depend on 



the initial state, and that do not require optimizations over all 
possible Kraus representations. Indeed, convenient classes of 
Kraus operators may be chosen inspired by physical consid- 
erations regarding the process under investigation. The power 
of this method was exemplified within the framework of two 
important problems in quantum metrology, the estimation of 
phase in optical interferometry and of transition frequencies in 
atomic spectroscopy. In optical interferometry, it captures the 
main features of the transition from the Heisenberg limit to 
the asymptotic shot-noise-like behavior. Furthermore, it im- 
poses a severe restriction on the asymptotic behavior of the 
estimation error: even for weak noise, the improvement on 
the shot-noise limit is at most by a constant factor, which does 
not change the 1 / ^/N behavior. 

Due to the omnipresence of non-unitary dynamics and the 
problem of parameter estimation in several fields of science, 
we envisage that this approach might find useful applications 
to other kinds of systems, involving for instance stochastic 
processes 113 ill and dynamical evolutions that depend non- 
linearly on the number of resources. 
Methods 

Outline of the derivation of the upper bound to the quan- 
tum Fisher information 

Given an initial pure state po = of a system S and 

an arbitrary process, which changes the state to p(x) = 
J2e ^h(x)po^-l{x), we expand the Hilbert space, introducing 
an environment E so that the total state in S + E undergoes a 
unitary evolution, described by 

\V{x)) = UsA*Ws\Q)b = Y,M*Ws\I)e. (15) 



Here \Q)e is the initial state of the environment. The states 
\£) e form an orthogonal basis in E, which is independent of 
the parameter x to be estimated and the Kraus operators 11^ (x) 
act on S. 

The upper bound for the quantum Fisher information is ob- 
tained from the inequality 



Tq = max F ( E 

£< S) «1 <E) 



(£f>®l< B >) 



< max F E n 

MS.E) \ J 



Cc 



(16) 

where F is the Fisher information defined in Eq. ([T). This 
inequality results from the fact that, for Tq, the maximization 
is made for all POVMs contained in the S space, while for 
Cq the maximization is for all POVMs in S + E. Therefore, 
Cq is an upper bound for Tq. The right-hand side of ( TTol l can 
be explicitly evaluated, since one is dealing in this case with a 
unitary evolution. It is then straightforward to obtain (O. 
Saturation of the upper bound to the quantum Fisher in- 
formation 

We show now that it is always possible to choose a set of 
Kraus operators such that the inequality in Eq. ( fT6l ) is trans- 
formed into an equality. The proof is based on Uhlmann's 
theorem 13211 and on the relation between Bures' fidelity 

Fb(pi,P2) = Trv/ p\' 2 p2p\ and the quantum Fisher in- 



6 



formation lfl2tl . Uhlmann's theorem implies that fioll 

(F B [p(x lca \),p(x)}) 2 = max|^ (:r) )|($(a; roa i)|*(x))| 2 , 

(17) 

where |$(a; rea i)) is an arbitrary purification of p(x lca \) in S" + 
E and the maximization runs over all purifications of 
p{x), also in S + E. On the other hand, Ref. [ 12] shows that, 
up to second order in 5x — x — a; rca i, 

(F B [p(x rcal ),p(x)}f = 1 - (Sx/2) 2 Tq [p(Zreal)] ■ (18) 

Expanding the right-hand side of WIS up to second order in 
5x, and comparing the resulting expression with the right- 
hand side of d 1 8t >. one gets (see Supplementary Material) 
(O. Therefore Cq coincides with the quantum Fisher infor- 
mation for some choice of Rraus operators. This proves the 
attainability of our bound. We show, in the supplementary 
material, that there is, in fact, an infinite number of Rraus rep- 
resentations that satisfy the above relation. 
Upper bound for quantum Fisher information in lossy op- 
tical interferometry 

A convenient choice for the Rraus operators corresponding to 
an interferometer with losses in one of the arms is 



tl e (9;a) = yii^e^-^W, (19) 

where a is the annihilation operator corresponding to the 
dispersive-arm mode, and r/ quantifies the photon losses (from 
i] = 1, lossless case, to rj = 0, complete absorption). The pa- 
rameter a defines a family of Rraus operators and it is used to 
minimize the value of Cq . This choice of Rraus operators is 
inspired by physical considerations. For a = (a = — 1) the 
process can be interpreted as a probabilistic photon absorption 
event, simulated by a beam splitter with transmissivity 77 and 
after (before) it a phase shift 9. If the beam splitter is placed 
after the dispersive element, then monitoring the environment 
leads to full recovery of the information, yielding therefore 
a non-reliable bound, which coincides with the one for the 
lossless case, while only partial information is obtained if the 
beam splitter is placed before the dispersive element (since 
in this case the deflected photons do not carry phase infor- 
mation). This last situation corresponds to the experimental 
setup used in l25ll . It may lead to a bound better than the one 



corresponding to the first choice, although not as good as the 
one obtained by optimizing over all values of a. 

The Rraus operators given by Eq. ( fT9] > lead to (see Supple- 
mentary Material) 

Hi(6;a) = 77 (1 - 77) (1 + a) 2 7i 

+ [l-(l-77)(l + a)] 2 n 2 , (20) 
# 2 (0;a) = [l-(l-77)(l + a)]n. (21) 

We get then, for all real values of a, 

C Q (p ;a) = 4[l-(l-ry)(l + a)] 2 (An 2 ) 

+ 4,7(1- 77) (l + a) 2 (n) , (22) 
where (n)o is the initial average number of photons in the 
dispersive arm, and (An 2 )o is the corresponding variance. 

When a = —1, one gets the trivial result Tq < 4(An 2 )o, 
while for a = 0, Tq < 4ry [r](An 2 ) Q + (1 - 77) (n) ] . The 
factor 77 before the brackets in this last expression comes from 
the energy loss, a classical effect. This expression implies 
that, only super-poissonian states ((An 2 )o > (n)o) may lead 
to an improvement of the phase estimation over the standard 
quantum limit, which corresponds to Tq — Arj(n)Q. 

For each value of (Att, 2 )o and (tt,)o, the minimization of Cq 
with respect to a yields ( fTTT i - see Supplementary Material. 

One should note that there is a choice of Rraus operators 
for which H2 = in (1201 1. and at the same time the term 
proportional to h 2 in Hi vanishes: the one that corresponds 
to taking in < f2Qb 1 + a = 1/(1 — 77). This implies, according 
to the discussion in the Supplementary Material, that, even 
in the presence of feedback 133143511 . the scaling in the phase 
uncertainty cannot be better than the standard quantum limit. 
Even though Refs. i33l435ll deal with global estimation, 
our method still applies, since (O is a lower bound for the 
uncertainty in global estimations, which is given by the 
Holevo variance QSQ . 
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